Probing the spin states of three interacting electrons in quantum dots 
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We observe a low-lying sharp spin mode of three interacting electrons in an array of nanofabricated 
AlGaAs/GaAs quantum dots by means of resonant inelastic light scattering. The finding is enabled 
by a suppression of the inhomogeneous contribution to the excitation spectra obtained by reducing 
the number of optically-probed quantum dots. Supported by configuration-interaction calculations 
we argue that the observed spin mode offers a direct probe of Stoner ferromagnetism in the simplest 
case of three interacting spin one-half fermions. 



PACS numbers: 73.21.La, 73.43.Lp, 73.20.Mf, 31.15.ac 

Systems of three charged particles interacting by 
Coulomb forces are the building blocks of a large variety 
of correlated quantum phases. An example of current 
topical interest is represented by the unusual fractional 
quantum Hall states with non-abelian excitations!^ On 
more general grounds, the interaction among at least 
three particles is essential to account for Stoner ferro- 
magnetism (SF), the paradigm of the tendency of itin- 
erant electrons, like those responsible for conduction in 
metals, to align their spins at the expense of their kinetic 
energy. 3 

Indeed, the simplest finite-size version of SF display- 
ing the crossover between normal (unpolarized) and spin- 
polarized degenerate (ferromagnetic) ground states at a 
threshold interaction strength is offered by three inter- 
acting fermions of spin one-half in two dimensions, as: 
(i) in one dimension the ground state is never spin polar- 
ized for any interactional (ii) the two-body ground state 
is always a spin singlet W 

Cold Fermi atoms confined in optical traps have been 
proposed as finite-size simulators of three-body inter- 
action phenomena 6 and SF physics. They are simple 
systems able to drive SF through the tunable short- 
range atom-atom interaction.^^ However, this capability 
is limited by the losses due to three-body recombination. 

Here we realize and study a solid-state version of the 
quantum simulator of three interacting fermions based on 
electrons confined in a semiconductor quantum dot (QD). 
We show that unpolarized and ferromagnetic states can 
be probed by monitoring the low-lying neutral spin and 
charge excitations. We recall that such collective modes 
in semiconductor QDs can be accessed by resonant in- 
elastic light scattering JEBHH We argue that QD systems 
with three electrons constitute a versatile emulator of the 
physics of SF. 

So far studies on low-lying neutral spin modes in QDs 
were performed by optical methods on nanofabricated 
AlGaAs/GaAs QDs and in self-assembled InAs QDs. In 
these experiments, however, ensembles of many QDs were 
investigated owing to the very low signal-to- noise ratio. 



For example, arrays composed of 10 4 nanofabricated Al- 
GaAs/GaAs QD replicas were studie d by us by means 
of resonant inelastic light scattering jnHl31i5] j n these sys- 
tems, however, inhomogeneities in the electron number 
distribution among the QDs contribute to a significant 
broadening of the excitation peaks in the detected spec- 
tra, leading to a systematic uncertainty in the identi- 
fication of specific contributions arising from different 
electron populations. Indeed, the emulation of SF re- 
quires access to sharp collective excitations of three elec- 
trons that are sufficiently isolated from excitations corre- 
sponding to other electron number configurations. This 
is achieved here in arrays of nanofabricated QDs which 
are diluted enough to suppress inhomogeneities linked to 
number fluctuations. We find that when the number of 
QD replicas probed in the light scattering experiments 
is decreased from 10 4 to 10 3 , the inhomogeneities are 
largely suppressed, leading to a single and sharp spin 
mode in the excitation spectrum. 

We argue in the following that this low-lying mode 
corresponds to a spin excitation of the three interacting 
electrons from the ground state with M=l, 5 = 1/2 
(being M the total orbital angular momentum and 5 the 
total spin of the system) to the excited state with M = 1, 
5 = 3/2. The assignment is corroborated by accurate f ull 
configuration interaction (CI) calculations! 11 * 13 * 15 * 17 ^ ^ 
The latter enable to link the measured excitation to the 
spin-polarized state with M = 0, 5 = 3/2 predicted by 
SF, which is an excited state at the actual value of the 
electron density. 

The three-electron system exhibits a transition to a 
ferromagnetic state as a function of the Wigner-Seitz pa- 
rameter r s , which is the radius of the circle whose area 
is the average area per electron, in units of the Bohr ra- 
dius. By knowing both the energy of the transition and 
the value of the energy spacing Jtljq among the shells of 
the potential trap, as retrieved from the measurements, 
we can thus deduce the energy difference between the 
absolute ground state of the system and the fully spin- 
polarized ground state. As it is shown in Fig. [I] this 
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makes it possible to place our results (filled circles in 
Fig. [lj on a phase diagram showing the path to the real- 
ization of SF. 

We recall that the essence of SF is summarized by the 
'Stoner criterion'.^ According to that, the ground state 
is ferromagnetic when n(Ep)J > 1, with u(Ef) being 
the density of states resolved at the Fermi energy and 
J being proportional to the exchange field that splits 
the energies of electrons of opposite spin. Hence two 
parameters, namely J and n(-Ep), control the relevant 
physics of the bulk. 

The finite-size analog of the criterion, derived in a sim- 
ple Hartree-Fock picture, is J/(Hujo) > 1, with l/(hu)o) 
replacing the density of states n(Ep). In our quantum 
simulator, the first parameter, Hujo, is obtained through 
the joint measure of both the electron number N in the 
QD and the pristine electron density n of the quantum 
well from which QD arrays are nanofabricatedP^ For 
fixed N, hojQ is decreased (r 3 is increased) by decreasing 
n. The last parameter, J, is measured from the lowest 
spin excitation, that is 2Hujq — J (cf. inset of Fig. [3]). 




FIG. 1. (Color online) CI energies of the unpolarized (red 
[gray] circles) and spin-polarized (black circles) three-electron 
ground states vs r s . The filled circles point to the values in- 
ferred from the inelastic light scattering measurement. Inset: 
zoom around the crossover. Ten harmonic-oscillator shells 
were considered in the CI calculation. A square QW of width 
25 nm and height 250 meV confines the motion in the out-of- 
plane direction. 



scattering measurements were performed at T = 2 K with 
a tunable ring-etalon Ti:Sa laser impinging with a spot 
size of 0.1 mm at normal incidence on the sample. The 
scattered signal from the QDs was collected by a triple 
grating spectrometer coupled to liquid N 2 -cooled multi- 
channel CCD. From previous experiments, it follows that 
a number of electrons ranging from 2 to 6 are expected 
to be confined in such QDs with a confinement energy of 
hu « 4 meV.™m 

We focus first on the sample with 10000 QD replicas. 
Figures [2^a) and [2^b) show the results of light scatter- 
ing measurements of spin (SDE) and charge density ex- 
citations (CDE), respectively, obtained with depolarized 
(perpendicular incident and scattered laser polarizations 
to detect the spin signal) and polarized (parallel incident 
and scattered laser polarizations to detect the charge sig- 
nal) configurations. Two prominent peaks emerge in the 
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FIG. 2. Top panels: experimental spectra of spin (a) and 
charge (b) density excitations detected at T = 2 K by means 
of depolarized and polarized resonant inelastic light scatter- 
ing, respectively. Bottom panels: Corresponding CI spectra 
obtained for N = 3 (solid line) and N = 4 (dashed line). 



A series of QD arrays were fabricated on a 25 nm mod- 
ulation doped GaAs/Alo.iGao.gAs quantum well (QW) 
by means of state-of-the-art electron beam lithography 
and inductively coupled-plasma reactive ion etching. The 
electron density of the two-dimensional electron gas was 
n = 1.1-10 11 cm- 2 , and the mobility 2.7-10 6 cm 2 /Vs. Pil- 
lars with diameters of 320 nm and aspect ratio of ~2/5 
were produced in 0.1 mm-sized square arrays, in which 
the number of QD replicas was varied from 10000 down to 
~ 1000 to explore the impact of inhomogeneous broaden- 
ing on the collective excitations. Resonant inelastic light 



spin channel at energies of about 4.4 meV and 5.6 meV, 
respectively, with the same intensity and full width at 
half maximum (FWHM) of approximately 1 meV. This 
value matches with the results of previous experiments, 
in which an a vera ge population of four electrons per dot 
was identified! 11 1 13 1 

In order to single out the contributions linked to spe- 
cific electron populations, we applied the full CI approach 
to retrieve the theoretical spectra linked to a series of 
electron populations, ranging from N = 2 to N = 6, 
with hujQ being fixed by the experimental value of n. This 
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analysis, which has spanned the wider expected range of 
confined electron numbers, enabled the identification of 
two prominent contributions to the experimental spectra 
arising from N — 3 and N = 4 electrons. Figures [2|c) 
and |2|d) show the agreement with the experimental data 
obtained in the case of three and four electrons per QD. 
In both spin and charge channels, the interplay of the 
contributions from these two electron configurations re- 
flects in a wider peak at higher energy. In the SDE spec- 
trum, in particular, there is a strong feature associated to 
the three-electron case at around 5.5 meV, which overlaps 
with the contribution of the four-electron configuration. 
Higher-energy peaks resolved in the theoretical simula- 
tion are not visible in the experimental spectra proba- 
bly due to their weaker intensities. The evolution of the 
SDE experimental spectra as the number of QDs in the 
illuminated array is progressively decreased is shown in 
Fig. [3] These data confirm the interplay between these 
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FIG. 3. (Color online) Resonant inelastic light scattering 
spectra of spin excitations obtained at T — 2 K for arrays 
containing a varying number of QDs. From top to bottom, 
the number of QDs in the array decreases from 10000 to 1110 
(the inter-dot distance increases from 1 /xm to 3 //m). Gaus- 
sian fits are superimposed to the experimental spectra in cor- 
respondence of three- and four-electron contributions (blue 
[dark grey] and red [gray] curves, respectively). The resulting 
global fit is shown as a green [light gray] curve on each spec- 
trum. On the right hand side, a scanning electron microscope 
image of the studied QDs is shown, in correspondence to each 
case. 



two different electron populations allowing to isolate the 
specific contribution of the QDs with three interacting 
electrons. The measurements carried out with varying 
the density of QDs allow us to confirm the broadening 
of the high-energy peak as mainly due to the simulta- 
neous presence of contributions associated to a different 
number of electrons, namely TV = 3 and N = 4 as mod- 
eled by the theoretical analysis. In Fig. [3J Gaussian fits 
centered at Raman energies matching with the CI predic- 
tions [Fig.[2jc)] have been considered in order to globally 
reproduce the experimental data. It follows that: (i) 
the two slightly shifted peaks relative to N = 3 (blue 
[dark gray] curve in Fig. [3| and N = 4 (red [gray] curve) 
present at ~ 5.6 meV broaden the overall parent struc- 
ture (green [light gray] curve) to ~ 1 meV FWHM for 
10000-QD array [Fig. ptta)]. (ii) the progressive narrow- 
ing of the peaks due to the vanishing of the N = 4 fea- 
ture in the less dense arrays makes the linewidth of the 
high-energy peak decrease down to the limiting value of 
0.4 meV FWHM for 1110 QDs [blue (dark gray) curve 
in Fig. [3jd)]. We thus eventually managed to realize a 
system in which only the population of three electrons is 
present. The spectrum of spin excitation measured for 
the array of 1110 QDs is reproduced in Fig. [4] together 
with the pertinent CI prediction for N = 3. 
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FIG. 4. (Color online) (a) Experimental resonant inelastic 
light scattering spectrum (black curve with points) and Gaus- 
sian fit (continuous blue [dark grey] curve) for the spin channel 
for the 1110-QD array. A scanning electron microscope image 
of the relative QD array is also shown, (b) Configuration- 
interaction calculations of spin density excitation for N — 3. 
The inset shows the energy diagram of relevant transitions. 
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By comparing the top and bottom panels of Fig. [4] 
and Fig. |3J it emerges that the N = 3 contribution has 
been identified and selected out of a more complex back- 
ground. The lowness of the signal, at the limit of the 
signal-to-noise detection threshold, implies that only the 
most intense peak, centered at ~ 5.5 meV, can be safely 
identified. 

The transition energies involved in the measurement 
are indicated in the inset to Fig. |4j The retrieved 
monopole SDE links the ground state, characterized by 
M = 1 and S — 1/2, to an excited state with M = 1, 
S = 3/2. The final state of this transition is also acces- 
sible via a Kohn mode excitation from the ferromagnetic 
ground state with M = 0, S = 3/2. Remarkably, the two 
spin-polarized states with M = and M = 1 differ only 
in the center-of-mass motion, whose excitation — exactly 
of energy Hujq — is insensitive to interactions. Therefore, 
with regards to the non-trivial part of the wave function 
that is affected by interactions, the final state of the ob- 
served SDE is a replica of the ferromagnetic ground state. 
This attribution is crucial to address the realization of a 
spin-polarized state relevant to the study of Stoner fer- 



romagnetism. 

In conclusion, we have identified the low-lying modes 
of three interacting electrons in a quantum dot by means 
of inelastic light scattering. By decreasing the number 
of optically-probed QDs we were able to suppress the 
inhomogeneities related to the simultaneous presence of 
different electron populations leading to a spin excita- 
tion spectrum dominated by a single and sharp (FWHM 
of 0.4 meV) peak. By a detailed comparison with calcu- 
lations based on a configuration interaction method we 
have linked to observed peak to a transition from the 
M = 1, S = 1/2 ground state to the M = 1, S = 3/2 
excited state of the three interacting electrons. These re- 
sults mark a starting point for the simulation of Stoner 
ferromagnetism with a three electron quantum dot and 
for the investigation of the competition between spin- 
unpolarized 5=1/2 and spin-polarized S = 3/2 states. 
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